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Abstract 
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We analyze the production of a single top squark at hadron colliders: pp/pp — ► t\ + X. 
The total cross sections and the transverse momentum distributions are presented in next-to- 
leading order QCD. The higher-order corrections render the predictions theoretically stable 
^ with respect to variations of the factorization and renormalization scales. The transverse 

momentum distribution is resummed to estimate effects of the small transverse momentum 
regime on possible analyses. Since the corrections increase the cross sections and reduce 
the theoretical uncertainty, the discovery range for these particles is extended in the refined 
\ analysis. 
00 ■ 

Hadron colliders like the Tevatron and the LHC will soon be able to either discover or strongly constrain 
physics at scales well above the Standard Model masses. From an aesthetic point of view, the most attractive 
realization of supersymmetry could be the minimal supersymmetric model (MSSM). However, R parity 
conservation is an ad hoc assumption, invoked to bypass problems with flavor-changing neutral currents, 
proton decay, atomic parity violation and other experimental constraints. From a more general point of 
view, these observables put tight limits on some, but not on all, R parity violating couplings. For collider 
searches, exact R parity predicts that supersymmetric particles can only be pair-produced. Limited by the 
beam energy, the Tevatron tends to run out of supersymmetry discovery reach if the strongly interacting 
squarks and gluinos become too heavy to be produced in pairs. But many models based on unification 
scenarios at some high scale prefer exactly this kind of mass spectra. If, in contrast, R parity is not an exact 
symmetry, the signals from single superpartner production can be nicely extracted in the low background 
environment at the Tevatron [[1|-|3]] or even at the LHC [Q] . 

In the case of a light scalar top squark, the baryon number violating coupling \'^ k induces the production 
process djd k — > t [Q]. It stems from the superpotential contribution |Q] 

W = X^Dpt, (1) 

where D c and U c denote charge conjugate right handed quark superfields. The couplings including at least 
one third generation flavor index {i, j, k} are currently only constrained in the combination | Ag^A^* 
The relevant Lagrangean readsQ 
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Cx>> = -2e a ^K jk 



tRad^Pndky + h.c. H (j < k) (2) 
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(a) a 0>K ^cx^ (c) 



Figure 1. Basic diagrams for the single stop production athadron colliders in antiquark-antiquark collisions: (a) 
generic diagrams of virtual QCD corrections; (b) generic diagrams for real gluon emission; (c) SUSY-QCD corrections 
to the stop propagator, including the non-Standard Model contributions to the running mixing angle. 

If we assume that the R parity violating coupling is a small parameter A" ^ 0.1, and the stop is not too 
light, the produced single top squark has a sizeable branching fraction into the 'classical' supersymmetric 
channel t — > bx + , where the light chargino decays to Ivxi ©• These decays have been calculated in 
next-to-leading order []8p ; as expected for a decay mediated by the weak coupling constant, the corrections 
is moderate (K < 1.2), and the scale dependence is well under control. Furthermore, the higher order 
corrections leave the dependence on the stop mixing angle essentially unchanged. 

Next-to-Leading Order Cross Section 

In leading order the partonic cross section for the production of a single light top squark ti with the mass m 
is given by 

olo = K qq - sin 9 t A | Oil = K qq — - o 1 ; K, 



qq w? l i - " M m 2 I a ' qq 4iVf ... 

\ / \ / c (3) 

where the stop mixing angle 9 t measures the fraction of the right handed squark field in the light mixed 
state. During the following analysis we will omit this over- all factor, i.e. we assume t\ = Ir. All cross 
sections can be scaled naively; we present the details of this argument below. 

The 0(X" 2 a s ) cross section includes virtual gluon contributions, as well as real gluon emission and the 
crossed channel with an incoming gluon (Fig. [j]). The infrared and ultraviolet divergences are extracted us- 
ing dimensional regularization. After adding virtual and real gluon contributions, and after renormalization, 
only collinear divergences remain. They are absorbed into the definition of the parton densities through 
mass factorization. Numerically we use the kinematical delta distribution in the Born cross section, eq®, 
to perform the convolution of splitting functions and the leading-order cross section. 

Diagrams involving heavy supersymmetric particles, like stops and gluinos, contribute to the virtual 
next-to-leading order matrix elements. From stop pair production [[Hip, we know that the contribution of 
these loops is numerically well below the remaining theoretical and experimental uncertainties. Therefore 
we disregard them in the following analysis to remain maximally independent of the underlying model. 
In the case of single stop production there are no vertex diagrams involving gluinos/heavy stops. The 
supersymmetric partners contribute through wave function renormalization of external quarks and stop legs. 
A second class of corrections arises from the Feynman diagrams shown in Fig. [j](c): the external light stop 
can mix to a heavy stop, which then couples to the incoming quarks. This is not part of the wave function 
renormalization, but it can be absorbed into a running stop mixing angle 8 t (fi), evaluated at the external 



2 The coupling ddt is in fact related to e.g. ddt by supersymmetry, but since the latter is a three scalar soft breaking 
parameter we can assume \'~^~ <C A^. 
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particle's mass scale. The Green's function with an external light stop leg will only depend on the heavy 
stop mass through the running mixing angle, i.e. the single stop production cross section will be proportional 
to sin 9 t in leading and in next-to-leading order. The numerical effect of the running mixing angle has been 
shown to be negligible [|]] . 

The hadronic pp and pp cross sections are obtained by folding the partonic cross section with parton 
luminosities. In leading order only qq initial states contribute to the production process. Especially for the 
production close to threshold, the gluon luminosity at the Tevatron is expected to be small. At the LHC, 
however, the qg — > t+jet cross section can be large. To estimate the gluonic contribution, we compute the 
total cross section for t+jet production cross section (Tab. 0). The transverse momentum is cut at a minimum 
value of 5 GeV. A resummed cross section, as it is derived later in this paper, cannot be used, since in the 
given order it does not distinguish clearly between the different incoming states. For a small stop mass of 
200 GeV, 40% of the t+jet events at the LHC involve an initial state gluon. This fraction drops to 25% for 
a stop mass of 500 GeV. Even at the Tevatron, light stops in association with a jet are produced through 
incoming gluons in 25% of the cases. This large fraction corresponds to the fact that for qq initial states the 
large valence quark luminosity does not contribute. 

The next-to-leading order total cross sections are presented in Figure ^. The leading order results are in 
agreement with Ref. [|2|], after taking into account the large effect from switching the parton densities from 
CTEQ4M to CTEQ5M1. The next-to-leading order cross sections are parameterized by the factor K = 
c"nlo / °"lo From the point of view of parton luminosities, both colliders are very similar; the initial state ds 
involves one valence and one sea quark, after summing over stop and anti-stop production. Therefore, the 
fraction of sea quarks and gluons in the proton dominates the K factor. For larger stop masses, the gluon 
luminosity at the Tevatron decreases rapidly; therefore the K factor drops to 1.15. At the LHC the change 
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Figure 2. Total cross section and K factors for the single tu production as a function of the mass. The couplings 
X" are all set to 0.1. The mixing angle dependence is omitted, sin 6 t = 1; the renormalization and factorization scales 



are fixed to the final state mass, and the parton densities used are CTEQ5L/M1 The curves are calculated for t 
plus t* production for the three different couplings X'Lf.. The dashed curves show the leading order, the solid curve 
the next-to-leading order results. 
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Figure 3. Renormalization and factorization scale dependence of the total cross section for a stop mass of 200 
GeV, normalized to the central scale value /x = m. Left: both scales are varied for two different initial states. The sb 
case looks very similar to the db curve. Right: factorization and renormalization scale are varied independently. The 
respective other scale is fixed to the final state mass. The set of LO and NLO curves are given for the Tevatron and 
the LHC. Note that the leading order renormalization scale dependence at both colliders is identical. All unspecified 
parameters are chosen as in Fig.Q 

is less dramatic. For the ds initial state the K factor starts from a slightly smaller value than at the Tevatron 
and slowly increases from 1.3 to 1.4 with increasing mass. 

One measure for the theoretical uncertainty coming from higher order corrections is the dependence of 
the cross sections on the renormalization and factorization scales. As a central value for both scales we use 
the final state mass (/i^ = iir = m). The dependence of the total cross section on the scales can be read 
off Fig. |3|: identifying both scales and varying them between m/2 and 2m shows how the uncertainty is 
reduced by the next-to-leading order prediction. In case of a ds initial state the uncertainty drops from 15% 
to 5% at the Tevatron and from 5% to even smaller values at the LHC. However, part of this unusually small 
scale dependence in leading and next-to-leading order is due to a cancellation between the factorization 
and the renormalization scale dependences. Especially at the LHC, the total scale dependence cancels out 
almost perfectly. Varying the scales independently yields an uncertainty of 10% in leading and 5% in next- 
to-leading order. 
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29.2 pb 


10.9 pb 



TABLE I. The total cross sections for single stop production alone, together with a gluon jet, and through the 
crossed gluonic initial state together with a quark jet. For the latter the perturbative cross section has been used, cut at 
a minimum value of 5 GeV. All parameters are set as in Fig. ||. 
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The Coupling A" in Next-to-Leading Order 

The R parity violating coupling ddtn is not an exclusive feature of the MSSM. The underlying model 
could as well be the Standard Model with an additional strongly interacting scalar; the coupling A" has 
no symmetry operation that connects it to any Standard Model coupling. The MSSM, however, adds new 
particles to the spectrum: the light-flavor squark, the left-handed scalar top quark, and the gluino can 
contribute to the next-to-leading order matrix element. For the process considered, the simple 'SM+t/j' 
model can be regarded as the low energy limit of the MSSM, after decoupling of the light-flavor squark, 
the heavy stop state and the gluino. In this case the diagrams involving any of these particles vanish; the 
projector on the right-handed stop state (sin 9 t ) becomes unity or zero, dependent on if the lighter stop is 
the partner of the left or the right-handed top quark. Since the A" couples strongly interacting particles, it 
has to be renormalized [0] . Starting out from the known renormalization constants with and without heavy 
supersymmetric particles 
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the relation between the four renormalization constants for the vertex leads to a renormalization and subse- 
quently a running of the R parity violating coupling A". Since the heavy particles decouple at typical hadron 



collider scales, the coefficient of A"(//fj) is determined by the low energy effective model: 



1 - 



3ocCp 1 



An 



e 



A"(Q 2 



1 _|_ ZosCp i ^ 



(5) 



The pole in the standard MS conventions is given as e/e = (47r) e r(l — e)/T(l — 2e). This treatment is 
analogous to the decoupling of heavy particles from the running of the strong coupling constant a s (fXji). 



The Cross Section for Small Transverse Momenta 

The perturbative calculation of the next-to-leading order cross section as presented in the previous sec- 
tion needs to be modified to investigate the px distribution of the top squark. If the transverse momentum 
becomes much smaller than the other scales in the process, i.e. the mass of the produced particle, the con- 
vergence of the perturbation series in a s deteriorates. The series in a s has to be replaced by a series in 
a s log 2 Q 2 /pt [0~0]- I* 1S possible to resum all terms at least as singular as p? 1 in the differential cross 
section da/dp T : the formalism has been successfully applied to Drell-Yan, single vector boson and vec- 



tor boson pair, heavy quark, and single Higgs boson production [|T5|-|i7|]. We start from the leading order 
asymptotic cross section in the limit of small transverse momentum: 



1 da^ ^ a s 
K qq — 



do dp T dy 



71 Spi 



171 
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=9,9 



+ I, <-> X. 



(6) 



Here x\ 2 = e ±y m/ y/S are the momentum fractions in the limit of small transverse momentum. The leading 
order coefficients for the single stop production can be extracted from the soft and collinear divergences in 
the perturbative result. In contrast to the heavy quark production case [|P7|], they do not consist of a Drell-Yan 
type contribution and additional final state radiation terms. Since the coupling A" violates baryon number 
there is no analogue without final state radiation, i.e. without strong interaction charge in the final state. The 
leading divergence in p? 1 is given by A = 2C F and B = —AC F . 
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The resummed differential cross section involves the integration over the impact parameter space of a 
Sudakov-like form factor: 



1 da r 
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(7) 



The integral in the exponent can be calculated analytically^; the boundary is canonically chosen as b = 
2e~ 7E . The Bessel function Jq{x) oscillates, and the amplitude numerically decreases only slowly for large 
transverse momenta. In this regime we use the expansion of the impact parameter integral for large px 
as compared to the lower limit of b in the Bessel period considered \$Wj\. The factorization scale for the 
resummation has been chosen as the mass of the final state particle, and the parton densities are evaluated at 
the point fi = bo/b. However, it has been shown that the resummed cross section is ill-defined for b > Aqcd; 
we follow one way of parameterizing nonperturbative physics by substituting the form factor W(b) by [JTJ3J] 



W(b) -> W(K) exp 



-b g 1 -b g 2 log 



2 i b max m 



1 + 



(8) 



The parameters g\,i can be fitted to Drell-Yan data |J15[], yielding g x = 0.15 GeV 2 and g 2 = 0.4 GeV 2 . The 
cutoff scale for the nonperturbative physics we choose as 6 max = (2 GeV) -1 . The numerical dependence on 
these parameters has been shown to be negligible [JT3|] . It is possible to match the large and small transverse 
momentum formulae using a matching function 
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(9) 



where the matching scale p™ atch ~ ?tt,/3 is the typical choice. For large transverse momentum the asymp- 
totic and the resummed cross sections give arbitrary unphysical results; numerically their contribution to 
the general cross section is faded out by the matching function. For small transverse momentum the per- 
turbative and the asymptotic function approach each other rapidly. The matching function smoothes out the 
behavior of the general solution in the transition region. The uncertainty induced by the particular choice of 
a matching function is only 0(a 2 ) [[15]]. 

The results of the resummation are presented in Figure |j. As expected from similar processes [JT3J] 
the resummed cross section peaks at pt ~ 5 GeV. The perturbative and the asymptotic cross sections 
numerically agree very well for small transverse momenta p T < 2 GeV. Qualitatively the perturbative and 
the resummed cross sections agree for intermediate p T > 10 GeV. However, quantitatively the matched and 
the perturbative cross sections can differ by up to 40% for p T < 60 GeV at the Tevatron and p T < 100 GeV 
at the LHC. As can be seen in Fig. |], the differential cross sections drop sharply with increasing transverse 
momentum; the choice of p™ atch should be taken into account as a theoretical source of uncertainty. 



3 This requires a truncation of the running of the strong coupling constant to leading order. We have fixed the leading 
order Aqcd to reproduce the correct value for a s (Mz); this choice introduces a numerical uncertainty, since the 
average renormalization scale in the process might be considerably larger, e.g. for a stop mass of 500 GeV. Moreover, 
the fitted value of Aqcd f° r the CTEQ5M1 parton densities assumes a next-to-leading order running. Changing the 
QCD scale in the given limits yields a numerical variation of < 5% for the resummed cross section. 
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Figure 4. Differential transverse momentum distribution with respect to the top squark. The four different line 
styles denote the perturbative (dot dashed), asymptotic (dashed), resummed (dotted) and general (solid) cross sections, 
as described in the text. The two sets of curves correspond to stop masses of 200 and 500 GeV. All other parameters 
are chosen as in FigQ 

Conclusions 

It has been shown [Q] that the search for single top squarks is a promising channel for the upgraded Tevatron. 
We have calculated the next-to-leading order signal and the t+jet production cross section. The total cross 
sections are enhanced by up to 30% at the Tevatron and 40% at the LHC. At both colliders, the crossed 
gluon-antiquark incoming states are important, since there is no pure valence quark production mechanism. 
However, the typical next-to-leading order features are similar to supersymmetric pair production processes. 
The factorization and renormalization scale dependence in next-to-leading order is reduced from typically 
10% to 5% or less. The variation of both scales simultaneously results in strong cancellations. For the 
exclusive t+jet production we have resummed the leading large logarithms for small transverse momentum. 
The matched distributions exhibit some numerical dependence on the matching; especially for large stop 
masses at the LHC the resummation significantly enhances the differential cross sections up to a matching 
point p™ atch ~ m/3. Finally, we note that if we had a way to predict the helicity structure of the coupling 
A", this process would currently be the only way to directly determine the stop mixing angle at a hadron 
collider [[nj. 
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